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Supposing the spin-triplet superconducting state of Sr2Ru04, the spin-orbit (SO) coupling 
associated with relative motion in Cooper pairs is calculated by extending the method for 
the dipole-dipole coupling given by Leggett in the superfluid ^He. It is shown that the 
SO coupling works only in the equal-spin pairing (ESP) state to make the pair angular 
momentum HL and the pair spin angular momentum id x d* parallel with each other. The 
SO coupling gives rise to the internal Josephson effect in a chiral ESP state as in superfluid 
A-phase of '^He with a help of an additional anisotropy arising from SO coupling of atomic 
origin which works to direct the d-vector into a5-plane. This resolves the problem of the 
anomalous relaxation of ^^0-NQR and the structure of d-vector in Sr2Ru04. 
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1. Introduction 

Nature of superconductivity of Sr2Ru04 has attracted much attention since its discovery 
by Maeno and his coworkers in 1994.^ Now it seems to have been accepted that the Cooper 
pair is in the spin-triplet state. ^ However, its gap structure of the triplet state has not been 
confirmed yet. Especially, the intrinsic direction of the d-vector is still under debate. 

The effect of spin-orbit (SO) interaction associated with relative motion of the Cooper pair, 
together with an atomic SO interaction and the Hund's rule coupling, is expected to be crucial 
for determining the stable superconducting state among possible states in the manifold of the 
triplet pairing states. It will be shown that effects of SO interaction are much more crucial 
than that of the magnetic dipole-dipole (DD) interaction which was essential for clarifying 
the nature of superfluid phase of liquid '^He.'^ 

A highlight at earlier stages of research of Sr2Ru04 is a result of NMR Knight shift by 
Ishida et aP which exhibits no decrease across the transition temperature Tc when the mag- 
netic field B of the order of Tesla is in the ab (RUO2) plane, suggesting that the Copper pairing 
is in the triplet state and the d-vector is in the direction parallel to the c-axis, perpendicular 
to the plane. This structure of d-vector was predicted by a group theoretical argument on 
the assumption that orbital and spin space are transformed together due to the "strong" SO 
coupling in forming the Cooper pairs. ^ However, the assumption of the "strong" SO coupling 
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is not self-evident. 

From experimental side, we should have been much more careful to draw a conclusion 
about the gap structure from this fact because the d-vector has a tendency to rotate in such 
a way that d and B are perpendicular with each other even within the aft-plane. Without the 
magnetic filed B or under the sufficiently low field, the direction of d-vector is determined 
by other small perturbation such as, DD or SO interaction (two-body and/or single-body), 
or the effect of sample boundary. The magnetic field of ~Tesla seems too large to draw the 
conclusion about an intrinsic nature of the gap, considering the low condensation energy of 
~ ksTc with Tc — 1.5K. Concerning this subtlety, we should remember the case of UPts, in 
which the Knight shift shows no decrease for all direction of the magnetic field B >0.5 Tesla,^ 
while it shows clear decrease across Tc when B <0.2 Tesla is applied along b- or c-direction, 
suggesting the intrinsic direction of d is in the foe-plane.^'''' 

Indeed, six years later, Murakawa et al reported that the Knight shift does not exhibit the 
decrease across Tc also in the case where the small field B (0.02 < B < 0.05 Tesla) is along 
the c-axis.^'^ This implies that the direction of d-vector identified by the former experiment 
is not intrinsic but is forced by the magnetic field. The most plausible interpretation of these 
Knight shift measurements made down to the low magnetic field of B 0.02Tesla is that the 
intrinsic direction of the d-vector is in the a6-plane and the anisotropy field in the a6-plane is 
smaller than 0.05Tesla at most. The latter conclusion is derived from the fact that the Knight 
shift does not exhibit any decrease across the Tc down to the magnetic field of B = 0.05Tcsla 
perpendicular to the c-axis.*^ Therefore, theories justifying the fact that d is parallel to the 
c-axis might lose their plausibility, and other explanations arc anticipated. Quite recently, 
it was shown that the intrinsic direction of d-vector can be in the a5-plane if the Coulomb 
interaction among electrons on the 2p orbitals of O (other than that on the 4d orbital of Ru) 
is taken into account^*^ together with the atomic SO interaction and the Hund's rule coupling 
at Ru site.^^ 

A role of the SO interaction between orbital and spin angular momentum of Cooper 
pairs may be crucial because it works to make the d-vector perpendicular to L, the 
pair angular momentum in the c-axis; therefore, the d-vector is in the a6-plane. This ef- 
fect may open a way to resolve another puzzle of anomalous NQR relaxation of ^^O in 
the superconducting regime. Indeed, the analysis shows that the dynamical spin sus- 
ceptibility J2^^''^Xzz{^.,|^)/|^\u)=uJtiQB exhibits a huge enhancement at < T < Tc while 
EqImXa;a;(q,w)/w|c^=c^^Qj^ and Eq ImX2/2/(q7 ^)/^U=u;nqr foUow the TMaw at < T < Tc, 
the canonical T-dependence for the anisotropic superconductivity with a gap with line(s) of 
zero.^^'^^ A possible explanation for this behavior is that the internal Josephson effect arises 
through the SO coupling and gives rise to the excess NQR relaxation other than the con- 
ventional relaxation in the superconducting state due to the quasiparticles contribution. The 
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purpose of this paper is to develop a theoretical framework to explain the effect of SO in- 
teraction of Cooper pairs in the equal-spin-pairing (ESP) state, and to resolve the puzzle of 
anomalous NQR relaxation observed in Sr2Ru04. 

2. Spin-Orbit Coupling of Cooper Pairs in Triplet State 

The SO interaction associated with the relative motion of quasiparticles is given as follows: 

= -f ^ E E i-^ • X [(2^ - m - m]] , (2-1) 

I j^i ''J 

where /xb is the Bohr magneton, mband band mass of electron, m* effective mass of quasiparti- 
cles, and g is defined as ^ = Hes/fJ-B, jJ-eS being the effective magnetic moment. The difference 
of prefactor of pi from that of pj is due to the so-called Thomas precession by the effect of the 
special relativity. The appearance of the factor mband/"i* in (2-1) is due to the vertex correc- 
tion for angular momentum, mband/'^*«) derived on an extended Ward-Pitaevskii identity as 
shown in Appendix,^^ and that for the spin density, The renormalization amplitude 

"a" in the denominators is cancelled by the weight of quasiparticles "a" . 

In this paper, we restrict our discussions within the triplet manifold of ESP. By the 
procedure similar to that described in Ref.^ for the dipole interaction, the interaction (2.1) 
leads to the SO free energy F^o, for the Cooper pairs in the chiral state with the pair angular 
momentum hL, as follows: 

F,o = -5so(id X d*) ■ L, (2.2) 

where the coefficient gso depends on the details of the dispersion of the quasiparticles and the 

pairing interaction. In the spherical model of three dimensions (3d), is given as 

riband ^ . 4^ - 1 

9so = 9d — X 4 _ , (2.3) 

m* 5^ 

where is the strength of the dipolar coupling in the "ESP" -superconducting state in 3d. In 
the cylindrical model or in two dimensions (2d), g^o is given as 

9so = 9d ^x8^2— • (2-4) 

Hereafter, we derive expressions, (2.2) and (2.3) or (2.4), starting with the interaction 
Hamiltonian (2.1) which is represented in the second quantization as follows: 

"^baiid I I A A ^ 



H,o = J J dridra— -— ^Vi(ri)V't(r2)a^^<5^5 



• "(ri - r2) X (-in) ((25 - l)Vi - 25V2))] M^2)M^i)- (2-5) 



By introducing the relative coordinate r = ri — r2, and the center of mass coordinate R = 
(ri -|-r2)/2, (2.5) is reduced to 



2 

A'b '^band 
-"so = Z — 

ri m* 



J J dRdr-^Vi(R + r/2)^p\{K - r/2)a«^5^5 
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V'5(R-r/2)V'^(R + r/2X2.6) 



rx i-ih) (^(45 - 1)V, - ^Vk^ 



The free energy due to the SO coupling is given by the expectation value of Hgo, (2.6). As 
in the case of dipole-dipole interaction,^ we rely on the following decoupling approximation: 

{^iiR + r/2)V':^(R - r/2)V'5(R - r/2)V^^(R + r/2)) 

~ (Vi(R + r/2)V':t (R _ r/2))(V'6(R - r/2)V'/3(R + r/2)) 

= (Vi(r/2)V't(-r/2))(^5(-r/2)^^(r/2)). (2.7) 

Then, since (2.7) is independent of R, 

= {Hso) = -^^^("^9 - 1)^ / dr^^afiS^S ■ F;^{r)[fx i-ih)Vr)]Fs^{v), (2.8) 
where V is the system volume and 

In terms of the conventional notation, 

F„^(r) = i{aa2)afi ■ F{r), (2.10) 

(2.8) is expressed as 

2 /"I 
Fso = -^^^^^(45-l)Fy dr-3(7^^,5^5((7fea2)^^(c7^(72)^5 

F^{r)[fx{-m)Vr)mir). (2.11) 

With the use of the identity, 
(2.11) is reduced to 

i^so = -^^^2ie,ft(4^ - 1)V J drlF,*(r)[r X (-in)V,)],F,(r). (2.13) 

Since the vector pairing amplitude is the eigen function of the relative angular momentum, 
the following relation holds: 

[fx {-ih)Vr)]iFe{r) = hUF^ir). (2.14) 
Then, with the use of the conventional definition of the d-vector as 

F(r) = F{r)d, (2.15) 

(2.11) is expressed as 

i^so = "X^^^^^ - l)2ieakdld,hLiV J drl|F(r)|2. (2.16) 
Thus, the free energy due to the SO coupling is given as the following form: 

i^so = -gsoiidx d*) ■ L, (2.17) 
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where the couphng constant gso is expressed as 

5so = /.|^(4^-l)2F|drl|F(r)P 
where the pair ampUtude F{r) is given by its k-representation as 

F{r) = ^e''''"F(k). 
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(2.18) 



(2.19) 



Exphcit form of F{r) depends on the type of pairing and dimensionahty of space. Let us 
first examine the ABM state in 3d case: Then the pairing amphtude F{r) is given as^ 



F{r) = ^' 



The k-integration of the part including kx is performed as 



dk 



47r r kpr (kpr)^ 



cos kpr sin kpr 



(2.20) 



(2.21) 



The integration of the part including ky is performed similarly, leading to the expression for 
the pairing amplitude F{r) as follows: 



F(r) 



2 r 



cos kpr sin kpr 



kpr {kpry 

Then, substituting (2.22), the r-integration in (2.18) is performed as 



(2.22) 



/ 



drl|F(r) 



dr 



cos kpr sin kpr 



r kpr (/cpr)^ 



where we used the formula of definite integral 

cos t sin t 

As a result, the coupling constant gso, (2.18), is expressed as follows: 

,2i 



/•°° dt 

Jo T 



1 

4 



5so = /x|^(45-l)27r^^y 
m* 



(2.23) 



(2.24) 



(2.25) 



This result should be compared with that for the dipole-dipole interaction given by Leggett 
for the ABM state:^ 



(2.26) 



Therefore, the relation (2.3) holds. With the use of this coupling g^, the free energy due to 
the dipole-dipole interaction is given as^ 



Fa = --ga\{d-L)f. 



In the case of ABM state in 2d, the pair amplitude F{r) = F{p) is given as 



(2.27) 



(2.28) 
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where p = (x, y) is 2d vector in xy-plane. The k-integration for the first term in (2.28) is 
performed as fohows: 

_ i dJpjkFp) 
kp dx 

= -i^JiikFp), (2.29) 

where Jn{z) is the Bessel function of the n-th order. By a similar calculation for the second 
term in (2.28), the pair amplitude (2.28) is given as 

F{p) = ^=^^J,{kFp) (2.30) 
P 

Then, substituting (2.30), the r-integration in (2.18) is performed as follows: 
/dr-L|nr)P ^ jm^jyi^^. 
= 2/dp^|F(p)|^ 



= 47r^2 r^[Ji(feFp)]^ 
Jo P 

= 27r*^ (2.31) 



where we have used the formula of definite integral 







j[Jiit)r = 2- (2-32) 



Therefore, the coupling constant gso, (2.18), is expressed as 

5so = Ml^(43-l)47r*V, (2.33) 
m* 

leading to the relation (2.4). 

In the case of 2d system Sr2Ru04, the free energy of dipole-dipole interaction is given as 
follows:^^ 



3c 



(2.34) 

where a(=3.87A) and c(=6.37A) are the lattice constant of the primitive cell of Sr2Ru04.^^ 
3. Non-Unitary State due to Spin-Orbit Coupling 

The SO coupling (2.17) induces the non-unitary component of d-vector in general. The 
deviation from the structure of unitary pairing is determined by the balance of energy gain due 
to (2.17) and the loss of condensation free energy -Pcond- Although it is not easy to compare 
both effects at arbitrary temperature < T < Tc, it becomes rather easy at T ~ Tc, in the 
so-called GL region. The -Fcond in the ESP state in the GL region, is given by the GL free 
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1 /r\-n\ r / T\ A^ + A? 7/-(^) K Al + Af" 

^ (3.35) 



F =l(^\ A T\ Af + A| 7C(3) K A| + A 
2 Vde^ V tJ 2 ^16 (7rfcBTe)2 2 



where (dn/de) is the density of states (DOS) near the Fermi level of quasiparticles in the 



normal state, k = I'^kl^/d^k l^kP)^; and A^^) is the tt (ii) component of the gap 



matrix. In the unitary state where A-|- = A_j_ = A, minimizing (3.35) with respect to A, 
is given as follows: 



1 /dn\ 8 1 . , ^ 



In the GL region, the coupling g^j^ is given as 

2 



3 



= (I) 7(1) ^<'*-^='' M^-^^Mf (l - I) . (3.37) 

Then, considering the case of Sr2Ru04, we use relation (2.4) in 2d and obtain the ratio of gso 



and l-^condl ^ follows: 



47r/iB :r [ln(l-14/3cec)] 1 - ;^ ■ (3-38) 



Let us parameterize the gap matrix as 

A = ^£^f-^-^ ' ]. (3.39) 
This is equivalent to represent the equilibrium value of d- vector as, 

1 Tj 

dox = —^=, doy = i—==, doz = 0. (3.40) 

Vl + r/2 \/l + r/2 

This d- vector in the equilibrium state is shown in Fig. 1. Substituting expression (3.39) into 
(3.35), after some standard calculations, we obtain the loss of condensation free energy AFcond 
as follows: 



AFcond = I ^cond I ^-y _^^2\2 ' (^'^^^ 



Fso = -9so^^- (3.42) 



The energy gain due to the SO coupling (2.17) is expressed in terms of (3.40) as follows: 

Here we have assumed that the pair angular momentum L is along the z(c)-axis since we are 
considering the case of Sr2Ru04. 

Therefore, the total free energy F{r]) = Fgo + AFcond as a function of rf is given as 

The degree of deviation rj from the unitary pairing is determined by the condition dF{rj)/dr] = 



7/18 



J. Phys. Soc. Jpn. 



Full Paper 



Lllz 



II X 



d^'lly 



y 



Fig. 1. Structure of d- vector. 



0, which is exphcitly expressed as 

(1 - rf) [-5.0(1 + rf) + 4|Fr^* |r/] = 0. (3.44) 

In the case where the deviation from the unitary pairing is small, 7/^ ^ 1 can be neglected 
compared to unity, so that condition (3.44) is reduced to a simple form 

9so 



In the case ^so > condition (3.44) is reduced 

r/ = 1, 

which implies that dp _L dg'; non-unitary state is formed. 
The prefactor of (1 — T/Tc) in eq.(3.38) is estimated as 



(3.45) 



(3.46) 



"2-band (45-1) 



/ dn 



- 



[ln(1.14/3cec)]^ ~ 4.8 x 10" 



L '^band 



(3.47) 



vnT ■ \ de y ■ " m 

where we have used the following relations, 5 — 1, (dn/de) = m* jcKf? ^ with c = 6.4A,^'^ and 
we have assumed 1.14/3cec — 20. Then, ratio (3.38) is estimated to be 



I punit I 
I cond I 



4.8 X 10 



-4 '^band 



-1 



(3.48) 



Therefore, except for narrow temperature region near T Tc, the condition rj^ <^ 1 holds so 
that the relation (3.45) is valid. Since the energy gain due to the SO coupling is multiplied by 
a small factor rj, (3.48), as in eq.(3.42), in order that the SO interaction predominates over 
the DD interaction, eq.(2.27), in the low-temperature region, we need another mechanism to 
make the d vector within the a6-plane, as will be discussed in the next section. 
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The temperature region where the effect of SO interaction predominates over the DD in- 
teraction is severely restricted near as (1 — T/Tc) < 10^"^, taking into account the relations 
eqs.(2.26), (2.33), (2.34), (3.42), (3.45), (3.48), and mband/^J- - 2.9.^'^ Therefore, it is techni- 
cally impossible to observe the non-unitary state by any probe considering the broadening of 
the Tc itself. 

4. Origin of Anisotropy of d- Vector - Anisotropy Fields 



The free energy giving anisotropy of d- vector due to the magnetic fields in ESP state is 
expressed as^ 

AF^agn = ^X^iV^P^(d • H)^ (4.49) 

where Fq is the Fermi liquid parameter for the spin susceptibility; Xz = (l + -^o )"^/^b('^^/*^^)- 
In the GL region, this is reduced to 

,1-- 

1 T 2 

AFmagn = ~ i i JL XzH 



k1 + F^ 



Tc ,.2 f dn\ ^2 



fii(-r] H\ (4.50) 



K (1 + ^0^)2^^ \de) 
This should be compared with the SO coupling constant 

"^bandc(45 - 1) TT 2 ( 8 1 / , ^ ^2 r, n 1 /I fl / 1 ^ 1 ( A r.^^ 

where we have used relations (2.4) and (3.37). Therefore, 

1 "Iband,, , ^aN2 8(43 - l)7r3 /dn\ 2 ri n 1 /I « m2 



AF„^ - «^ (' + [Te ) ['"(l-l^A'o))' . (4.52) 

The actual competition between the effects of the SO coupling and the magnetic fields is given 
by 5so»?/AFmagn: In the region r/^ < 1, 

9soV 1 n , ™^2 27r2 / "^band (4^ - 1) . 2 /dn' 



AF 

'-^^ magn 



= 3.0xlO'^(l + F„«f (^)'r.= (l-|)"' [gauss-^l, (4.54) 

where we have used the estimation (3.47), 5 2± 1, (dn/de) = m*/cTrh^, with c = 6.37A,^^ and 
we have assumed l.lAPc^c ~ 20. With the use of experimental values of Sr2Ru04, T^. 2± 1.5K, 
rrib^nd/m ~ 2.9,^'^ and (1 + F§) 2± 0.5,^^ 

^g- ^ 1.4 X 10' ± (1 - I) [gauss-2]. (4.55) 
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Then, the anisotn 

5so^/AFmagn = 1: 



Then, the anisotropy field i?a°^^^ due to the two-body SO effect is given by the condition 



/ J. X -1/2 

iJf (2) ^ 12 X 10 f 1 - - J [gauss]. (4.56) 

In the hmit of T — Tc, tj approaches 1 so that the ratio (?so^/A-Pmagn is estimated as 

^^""^ ~ 0.25 X 106^2(1 + i^o")'^^^? [gauss-2] (4.57) 



~ 4.1 X 10^ ^9 [gauss-2]. (4.58) 
H 

Then, H^^"^^ is given as 

iff ~ 6.4 X 10^ [gauss]. (4.59) 

The dipole-dipole interaction also gives rise to the anisotropy of d- vector, as discussed 
in superfluid ^He^' The favorite direction of d- vector is along L, relative angular momentum 
of Cooper pairs, i.e., c-axis, in contrast to the effect of the SO interaction which works to put 
the d- vector in the ab-plane. By the analysis similar to the above, using expressions (2.34), 
(3.37), and (4.50), 

^^^^^^ 0.39 xl0^i^(l + Fo«)^^T,^:. 3.6x10^^ [gauss"^], (4.60) 

'-^-f^magn n Vfl £1 

where, in deriving the second relation, we have used c/a = 1.65, — 1.5K, m* /m ~ 16.0,^^ 
and (1 + Fq) ~ 0.5.^® Then, the anisotropy filed H^^ is estimated as 

H^^ ~ 1.9 X 10^ [gauss]. (4.61) 

This cannot predominate over H^^'^\ eq.(4.59), in a very restricted temperature region near 
Tc, as discussed in the last paragraph of §3. However, in the wide range of temperature 
< T < Tc, iff^ predominates over H^^'^\ Therefore, the favorite direction of d- vector is 
the c-axis as far as the SO interaction and the DD interaction of Cooper pairs (both of which 
are two-body effect) are taken into account, which is in contradiction to the Knight shift 
measurements,^'^ as discussed in the fourth paragraph of §1. 

Another mechanism giving the anisotropy of the d-vetor is an atomic SO interaction 
(single-body effect) together with the Hund's rule coupling.^^"^^ In particular, Yanase and 
Ogata showed,^^ on the basis of the 3rd order perturbation calculation of multiband Hubbard 
model for the pairing interaction,^^ that the favorite direction of d-vector is the c-axis and 
the anisotropy field H^^^^ due to the single-body SO interaction is of the order of 10~^Tc ~ 
2 X 10^ gauss. Then, adding H^'^, the total anisotropy field = H^'^ + i^f^^^ amounts 

to i^totai ^ 4 X 10^ gauss. This is also in contradiction to the Knight shift measurements,^'^ 
which shows that d-vector is perpendicular to the c-axis down to H = 200 gauss. 

Quite recently, however, it was shown by Yoshioka and the present author-*^^ that the 
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stable direction of d-vector changes from the c-axis to the a6-plane if we take into account 
the Coulomb repulsion of electrons on 2p-orbitals at O site"*^" which cannot be neglected as 
shown in a band structure calculation.^"^ Indeed, the anisotropy field H^^^^ is estimated as 
of the order of lO^^Tc by the calculations similar to ref.22. There is a possibility that this 
anisotropy field wins H^'^, ensuring that the favorite direction of d- vector is in the a&-plane. 
Furthermore, anisotropy of d-vector is in the a6-plane arises through the process breaking the 
conservation of z-coniponent of Cooper pairs by the atomic SO interaction.^^ It is noted that 
this process gives rise to a weak non-unitary component in the Cooper pair state. 

5. Internal Josephson Oscillations 



It turns out that the so-called internal Josephson effect due to the SO coupling is possible. 
When the d-vector in the equilibrium is given by eq.(3.40), its real and imaginary parts, dp 
and dp, are 

1 , .„ V 



dn = — , =x, dn 



(5.62) 



where x and y are the basis vector in the x- and y-directions. As the real component of 
d-vector, dg, deviates from equilibrium as shown in Fig. 1, the d-vector is 

cos 9 



d' 



Then, the diagonal components of the gap matrix are given as follows: 

1 



(5.63) 



-dx + idy 



(— cos 9 — 7] + isinO) 



V^rT2r7COsy+^ 



exp 



-itan 



sm( 



COS 9 + ri 



(5.64) 



and 



dx + idy 



: (cos 9 — T] + i sin 9) 



\/l — 2'q cos 9 + 7f 



exp 



itan 



-1 



sin0 



cos 9 — rj 



(5.65) 



Therefore, the phase difference between tt" ^'^d J^-component of gap matrix is calculated 
as 



tan 



tan 



-1 



sm ( 



cos 9 + T] 

sin 26* 

cos 26 — -rf 



+ tan" 



-1 



sm( 



cos 9 — -q 



Up to the 0{rj^), it is approximately given by 29: 

M = 29[l + 0{ri^)\. 



(5.66) 



(5.67) 
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For the smaU oscillations for which 9 <^ 1, the gap matrix is given up to 0{rf') as follows: 

\ (l + r,)e"J 

With the use of d- vector, eq.(5.63), the "pair-spin" is calculated as 

{id xd*)^ = 2((? X (F')^ = 2r/ cos 9. (5.69) 

Then, the free energy due to SO coupling (2.17) is 

Fso{9) = -gso{2ri)cos9, (5.70) 

where we have assumed that L = z as discussed above. 

The gap structure, eq.(5.68), and the free energy, eq.(5.69), have almost the same structure 
as those for the case of dipolar coupling in the A-phase of superfluid ^He.^ Only difference 
is that |A-|-|-| and |Aj_|,| are slightly different and has weak ^-dependence. Even though, the 
equation of motion for A^ and AN = N^ — N^, N-^(i) being the electron number with the spin 
t (I), are 

>t-iV,)^-i|^ = -^.in.. (5.71) 

and 

^(29) = 2A^ = ^(iVt - ^i)- (5-72) 

These coupled equations, eqs. (5.71) and (5.72), describe the harmonic oscillations whose 
angular eigen frequency is given by 

Xz 

9soil^B/hf gso (5 
Xz ^li'^condl 

where expression, eq.(3.45), has been used for 77 ^ 1. With the use of eqs. (3.38), (4.51), and 
= (1 + F^)~^iJ,^{dn/de), Cl'^, eq.(5.74), is expressed as 

Substituting eq.(3.47), we obtain the eigen frequency J7 as 



0~4.3x 10\/(l + Fo")/Krc-^^^ [sec-i]. (5.76) 



m 



6. NQR RelELxation Rate due to Internal Josephson Oscillations 



The dynamical uniform spin susceptibility in the A-phase in the ESP manifold has been 
established in the context of superfluid ^He by Lcggett and Takagi as:^^ 
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where the damping rate T is give as 

r = 7or02, (g_7g) 
where r is the hfetime of quasiparticles at normal state given as 

where 6 is a constant of 0{1). The coefficient 70 in eq.(6.78) is defined as 

70 = [1 - y(r)]-iy(T) — (e.so) 

where Y{T) is the Yosida function for the p-wave ESP state. ^ 
The NQR/NMR relaxation rate is given by 

I _ A ^ lmxz{q,uj) , . 

TYT-^l^ ^ ' ^ ^ ^ 

'^i^ q<qc 

where ^4 is a constant arising from the coupling between nuclear and electron spin fluctuations. 
Imaginary part of Xz{^)i eq.(6.77), is 

I"^X^(^) _ ^ Tor 



+ (7o<^t) 



2 



Since expression (6.82) is valid for the wave number smaller than inverse of the coherence 
length, the size of the Cooper pair, the cut-off wave number q* should be set as q* ~ r{Tr/^o), 
where r (<1) parameterize the cut-off size. 

Then, considering the case of Sr2Ru04 where w 2± 1 x 10^ and eq.(5.76), the NQR relax- 
ation rate due to the internal Josephson oscillations is given as 



where til is the 3d number density of lattice sites, and n2d is the areal number density 
of quasiparticles. Here, the factor (oj/Q)'^ has been neglected compared to unity, because 
Omega, eq.(5.76), is estimated as ~ 10^ ~ lOw using (1 + Fg) ~ 1/2,^^ k ~ 1, Tc ~ 1.5, and 
mband/'"^ — 2.9.^^ The ratio of ^0 and the lattice constant a in the plane is estimated as in 
the BCS model: 

^ = 1.1x10-1^. (6.84) 
With the use of the experimental values, ^0 = 1050A, and a = 3.87A,-'^^ Tp/Tc is given in turn 



as 



~ 2.5 X 10^. 



(6.85) 



13/18 



J. Phys. Soc. Jpn. Full Paper 



In the cylindrical or 2d model, the DOS is given as follows: 

'dn\ ^ 1 n2d 



de J c kBTp 

Then, the relaxation rate, (6.83), due to the internal Josephson effect is expressed as 



(6.86) 



M , loA^^^X.* V^V (6.87) 



This expression should be compared with the Korringa relation in the normal Fermi liquid 
state: 

q 

xo(q,w) 



^y^im 



q 



l + /oXo(q,u;) 



(6.89) 



EqImxo(q,a;)/a; 



AXz \dej Qc 2 2 

where x^(q,a;) is the transverse dynamical spin susceptibility in the normal state, Xo(q)^^) is 
the particle-hole propagator of quasiparticles (with xo(0) 0) = ^(dn/de), fcp is the Fermi wave 
number, and Qc 2± \/2^a is the wave number cut-off representing the lattice effect. 

In the case of Sr2Ru04, Xz — 2.0 x /x|(dn/de) or ~ —0.5,^^ 70 defined by eq.(6.80) is 
approximately given as 

70 ~ 2.o[i -y(r)]-iy(r). (6.92) 

The parameter lot in eqs. (6.82) and (6.87) is given in the present case, a; 2± 1.0 x 10^ sec~^, 
as 

UT ~ 7.6 X 10-5 X 6^ ~ 1.9 X 10-1^ , (6.93) 

where wc have used cqs. (6.79) and (6.85). 

The longitudinal relaxation rate of NQR, normalized by that in the normal state, eq.(6.91), 

is 

(i/rir)s(j, _ 6.5x10 ,,,,,, ^./liV 70 



*=4(1 + F„«)(|) (6.94) 



(1/TiT)n {qc/kY)a?n2d 
The Yosida function Y{T) in 70, eq.(6.92), is estimated in the standard manner by assuming 
that the pairing interaction Vk,k' = — Fcos (/?k cos (/^k being the azimuth in the k-space 
of a6-plane, and that the superconducting gap Ak = Acos(^k follows the weak-coupling 
gap equation. Since, in expression (6.94), there exist parameters h and r that are difficult 
to estimate microscopically, we choose them so as to reproduce the observed temperature 
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dependence of NQR relaxation rate.^^ Instead of b and r, two independent parameters can be 
chosen also as 

^ " ( ^il ^ 2° br', (6.95) 

and 

15 = 0.39;^. (6.96) 

In Fig. 2, we show the rcsTilts of the NQR relaxation rate in the superconducting state, 
(l/riT)g = (l/TiT)g(j) + (l/TiT)g(Q-), where (l/riT)g(Q^ is the quasiparticles contribution 
and is replaced by experimental values of (l/TiT)f„ for two sets of parameters, (I) C = 0.2, 
D = 0.44, and (II) C = 0.55, D = 1.0. Agreement of experimental measurements and the 
theoretical results, based on eq.(6.94), is rather nice, although the theoretical ones include ad- 
justable parameters and relatively crude approximations have been done. To our best knowl- 
edge, the unusual relaxation rate {l/TiT)c has not yet been explained. So, our theory may be 
the first one that explains the unusual behavior of NQR relaxation rate.^^ 

7. Summary 

We have obtained a formula for the spin-orbit coupling of the Cooper pairs in ESP state 
of Sr2Ru04, giving rise to the internal Josephson oscillations of d-vector in the a6-plane if 
the stable direction of d-vector is in the o6-plane. The latter condition is confirmed by a 
recent theoretical finding of Yoshioka and the present author^^ that the stable direction of 
d-vector is in the a6-plane on a realistic model of Sr2Ru04 by taking account of atomic spin- 
orbit interaction and the Hund's rule coupling among electrons on 4d orbitals at Ru site. The 
anomalous temperature dependence of NQR relaxation rate {l/TiT)c was explained by the 
theoretical formula, eq.(6.94), due to the internal Josephson oscillations of d-vector in the 
a6-plane that induces the oscillations of spin polarization in the direction of c-axis. 
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Fig. 2. NQR relaxation rate: Comparison of experiment and theory. Filled (open) circles represent 
data of measurements of NQR relaxation rate due to spin fluctuations along the c-axis (6-axis).^^ 
Dashed line represents (l/riT)s(Q) due to the quasiparticles contribution which follows the T^-law 
in a wide region T < T^. Solid lines (I) and (II) represent (l/riT)s = (l/riT)s(j) + (l/Tir)s(Q) 
for the parameter set of eqs.(6.95) and (6.96), (I) C = 0.2, D = 0.44, and (II) C = 0.55, D = 1.0, 
respectively. 



Appendix 

In this appendix we show how the factor mband/"*-* appears in (2.1) on the basis of 
an extended Ward-Pitaevskii identity. Suppose the system is subject to the low rotation 
5i}{r) which is slowly varying with respect to r. Then the term — (r x p) • is added to the 
Hamiltonian. Then, the variation of the Hamiltonian is given by the term 

- J dr^Ur) ([r X p] • 5Jl(r)) ^,(r) = - J drVi(r) ([dQir) x r] • p]) V'„(r) (7.97) 

where p = —iKV. By this perturbation, we obtain for the variation of the Green function: 

5G = -G(p)(iVp X p) • 5nG{p + k) + ^G{p)G{p + k) 

^r,^,a/3(p, q; k)G{q){i^, x q) • 6nG{q + k), (7.98) 
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where k = (k, 0) (we assume k to be extremely small). On the other hand, the addition of 
(7.97) to the Hamiltonian leads to the transformation of the momentum in the Green function 
as 

p ^ p - mband'^^^ X (iVp). (7.99) 

Hence, 



— = -mband(iVp) X — , (7.100) 



Since the relation 



SG ,^,dG 
^ = -mba„d(iV,) X — , 

as k ^ 0. Consequently, in the limit of dCl 0, and k ^ 0, we obtain for G{p) describing the 
quasiparticles: 

QQ—^ i f d^o 

mband(iVp) X ^— = -(iVp X P) + 2 y (^r^/3,a/3(P'9){G'(9)(iVg X q)G(g)}fe. (7.101) 

« * ' (7-102) 

op a m*a 

holds for the quasiparticles near the Fermi level, the relation (7.101) near the Fermi level is 
rephrased as 

-^(iVp X p) = -(iV^ X P) + ^ / ^^Ua0iP^l){G{qm, X q)G(g)},. (7.103) 

This explains why the vertex correction of spin-orbit coupling is given by mband/'^i*) leading 
to expression (2.1) after the factor 1/a has been cancelled with the renormalization amplitude 
a of quasiparticles. 
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